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TECHNICAL MEMORANDUM X-53996 


THE GENERATION OF X RAYS IN PLANETARY NEBULAE 

SUMMARY 


The evolution of planetary nebulae and the production of X rays will be 
summarized from the literature. Planetary nebulae emit continuous radiation 
which oi’iginates by free -free transitions involving kinetic energy losses of 
free electrons in the electrostatic fields of ions. These free-free emissions 
are responsible for the radiation observed in the radio frequency spectra. They 
also contribute in the visible and infrared regions. Therefore, if X rays are 
generated by planetary nebulae, it can be assumed that they will also be produced 
by free -free transitions. However, free -free radiation will be measured along 
with free -bound radiation. A formula is given for the intensity of the combined 
continuum for an ionized gas with a Maxwellian distribution of the electron 
velocities. 


INTRODUCTION 


In a statistical study which was undertaken to determine the degree of 
correlation existing between the positions of X-ray sources and classical 
cepheids, galactic novae, planetary nebulae, and Wolf-Rayet stars, the correla- 
tion coefficients obtained for the cepheids and Wolf-Rayet stars were small 
and those obtained for the galactic novae and planetary nebulae were high. It 
can be concluded that X-ray sources are more strongly associated with disk 
Population II objects than with extreme Population I objects. X rays are 
produced in supernovae as synchrotron radiation and in novae as deceleration 
(free-free) radiation. Planetary nebulae emit continuous radiation which 
originates by free-free transitions involving kinetic energy losses of free 
electrons in the electrostatic fields of ions. These free -free emissions are 
responsible for the radiation observed in the radio frequency spectra. They 
also contribute in the visible and infrared regions. Therefore, if X rays are 
generated by planetary nebulae, it can be assumed that they will be produced by 
free-free transitions also. Distant encounters between protons and electrons, 
involving small losses of energy, are more numerous than close approaches 
which produce large energy changes. Therefore, the free-free emission is 
more important at lower frequencies or longer wavelengths. 



THE EVOLUTION OF PLANETARY NEBULAE 


To maintain the energy balance in a planetary nebula, the star has to 
possess a degenerate core with a contracting outer region. The most highly 
evolved stars approach the white dwarf state. Since the mass lost during the 
formation of the planetary nebula occurs at an early stage, the surrounding 
shell is too faint to be detected in later stages. A significant number of the 
white dwarfs has evolved through the stage of planetary nebulae. 

To estimate chemical composition of a planetary nebula is: 

Element H He C N O F Ne Na S Cl Ar K Ca 
log N 12.00 11.25 8.7 8.5 9.0 5.2 8.2 6.0 8.0 6.5 6.9 5.8 6.2 

The planetary nebulae have an original mass of about 1. 2 with 

a lifetime, after nuclear reactions begin, of 7 x 10 9 years. The star-nebula 
systems observed as planetary nebulae represent the gravitational contraction 
phase following nuclear fuel burning in stars of M = 1.2 M^, Less massive 

stars do not have to lose mass in this phase and will not be identified as 
planetary nebulae. The more massive stars may evolve in such a manner that 
the mass loss occurs slowly, failing to produce the bright shells. 

An evolutionary track of log —■ against log T on the Hertzsprung- 
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Russell diagram is shown in Figure 1. The positions of radius R 0 = 0. 062 
and R = 0. 1, 0.2, and 0. 4 parsec are indicated. Since the nebulae are 
expanding, evolution takes place in the direction shown by the arrows. 

The electron density of a nebula, and the luminosity and temperature of 

its central star vary with the time t. Initially N is large, the luminosity 

© 

L equals 60 L , the temperature T equals 32 000° K, and the nebula is 

optically thick. The star evolves to L = 25 000 L , T = 60 000° K, the 

nebula expands, and the optical thickness decreases. At a time t 0 , radius 
R 0 , the nebula becomes optically thin. For t 0 < t, continues to decrease 

because of expansion. Then the temperature rises to 10 5 ° K at approximately 
constant L, and the optical thickness continues to decrease during this 
phase. Thereafter L decreases rapidly; the temperature remains at 10 5 ° K 
until the luminosity declines to 100 L . Finally, the decrease in luminosity 
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resulting from the onset of degeneracy becomes more important than the 

expansion, and the optical depth increases. At a radius E , the nebula 

i 

again becomes optically thick. E 0 = 0. 06 pc and R t = 0. 6 parsec. 

The results of Figure 1, the main sequence, the horizontal branch, and 
the white dwarfs are drawn in the H-R diagram of Figure 2. Using an expan- 
sion velocity of the nebula of 20 km sec -1 , a time of i. 7 x 10 4 years is obtained 
for the evoluation from R 0 = 0. 06 pc to R = 0.4 parsec. The entire track 
shown in Figure 2 is described in approximately 5 x 10 4 years. Tracks for 
two models, calculated by Hayashi and collaborators for Russell mixture stars 
evolving without nuclear processes, are also included. An initial contraction 
at constant luminosity is arrested by the effects of degeneracy and followed by 
cooling at constant radius. These models show some measure of agreement 
with the results for the central stars of planetaries, because the luminosity 
drop in the later stages of their evolution is a consequence of degeneracy, and 
because the final stage of their evolution points to the white dwarfs. The 
average mass of the central stars of planetaries, as well as of white dwarfs, 
is assumed to be 0.6 M Q . The average nebular mass is 0. 6 M . These stars 

belong to the disk population (Tables 1 and 2) . 

There exists no satisfactory explanation of the processes leading to 
mass ejection and to the increase of the radii of the central stars during the 
initial stage of their evolution. 


THE PRODUCTION 0FX RAYS 


The radiation emitted during encounters of electrons with ions in a 
fully ionized gas of low density is called "deceleration radiation. " The tei’m 
"free -free transitions" is customary for the same process in the astrophysical 
literature. The formulas for the deceleration radiation by a Maxwellian 
distribution of electrons will be given. 

The initial and final energies of a free electron may be represented by 
the expressions: 


E = ^mv 2 and E' = E - ha 
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The free-free emitted power per unit volume of a gas is obtained by 
integrating a quantity proportional to g/v over a distribution function propor- 
tional to the expression: 


e 


_E_ 

kT 


v 2 dv 


It is then found that the weighted average of g for a Maxwellian distri- 
bution may be expressed in two ways: 


g(T, r) = e 


hr E 

OO E 

kT r kT — 

J g(E, v) e d kT 

hr 


and 


kT — 

g(T, r) = g(E' + hr, v) e d kt 


The equation for the power density per sec -1 
emitted in free -free encounters is: 


/ 2tt m \ 2 


P (T, v) _ 8Z 2 e 6 

" 3(47re 0 ) 2 m 2 c 2 ^ 3kT 


4T 


N. N 
i e 


bandwidth per steradian 


g(T,r) e 


hr 

kT 


This formula yields the free -free emission. It contains the exponential 
factor, because free-free emission of frequency r can occur only if the initial 
electron energy exceeds hr. However, deceleration radiation can be emitted 
by electrons of any initial velocity. In free -bound radiation, the electron is 
captured by the ion. The two emissions will be measured together. 
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The free-bound radiation differs from the free-free emission by the 
quantization of the final energy. The inclusion of the free-bound emission will 
cancel the exponential factor in the power density formula, approximately. 

The successive entry of more free-bound continua yields a saw-tooth function 
which remains constant on the average, so that the exponential drop of the 
free-free emission is nearly cancelled by the free-bound contributions. At 
low frequencies or at high temperatures the free-bound contribution becomes 
negligible. Applying the relevant Gaunt factor, the emission spectrum from 
a gas for a combined continuimi is given by the expression: 


P(T, v) 
4tt 


= N. N 
1 e 


C v e 


hr 

kT 


g(T,n) + f(T, v) 


h 

kT 
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4t r Z 2 e 6 
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/ 8kT 
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X =m( v) 
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This formula gives the intensity of the combined continuum for an ionized 
gas with a Maxwellian distribution of the electron velocities. The values of 
(g + f) are given in Table 3 for 0 = 1, which is equivalent to T = 1. 58 x 10 5o K 
for hydrogen. The free-bound continuum becomes increasingly important 
from A = 1/i down. g(T, v ) , alone, is to be computed only if: 


hr << kT 


g(T, v) f(T, v ) is the quantity which has to be compared with experimental 
emission data. 

In the absence of computed values of f(T, v ) , g (T, v ) can be used if 
the exponential factor is dropped. 

o 

Evaluating the total power for 100 A, which is the long -wave limit for 
X rays, the following values are substituted into the equation: 


c 

V 

T 

N 

g + f 

joule m 2 

m sec 1 

°K 

m -3 


1. 82 948x1 0~ 45 

2.46936x10® 

1. 58x10 s 

4. 169x10® 

6. 58 

The result obtained is: 





P(T,n) 
47 r 


1.73 x 10 40 watt m 3 sterad 1 ( sec -1 bandwidth) 1 


6 




TABLE 1. OORT'S CLASSIFICATION OF STELLAR POPULATIONS [2], 


Population 

Halo 

Population II 

Intermediate 
Population II 

Disk Population 

Intermediate 
Population I 

Extreme 
Population I 

Some Typical 
Members 

Subdwarfs 
Glob. Clusters 
RR Lyr Var, 
With Per. > 0. d 3 

High-Vel. F-M 
Long-Per. Var. 

Planet. Neb 
Bright Red 
Giants 
Novae 

Weak-Line 

Stars 

Strong-Line 
Stars 
A Stars 
Me Dwarfs 

Gas 

Supergiants 
T Tauri Stars 

|z| (parsecs) 

2000 

700 

450 

300 

160 

120 

I Z | (km/sec) 

75 

25 

18 

15 

10 

8 

Axial Ratio 

2 

5 

=^25 



100 

Concentration 
Toward Center 

Strong 

Strong 

Strong 


Little 

Little 

Distribution 

Smooth 

Smooth 

Smooth 


Patchy, 
Spiral Arms 

Extremely 
Patchy, 
Spiral Arms 

Z h. 6 ' 

( Schwarzschild) 

0. 003 

0. 01 


0. 02 

0. 03 

0. 04 

Age ( 10 9 years) 

6 

6.0 to 5.0 

5 

1. 5 to 5 

0. 1 to 1. 5 

< 0. 1 

Total Mass 
(10 9 O) 

16 

47 

47 

47 

5 
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TABLE 3. AVERAGE GAUNT FACTORS FOR 0 = 1 [3] . 


log (v/v 0 Z 2 ) 

g + f 

log (v/v Q Z 2 ) 

g + f 

log ( v/v 0 Z 2 ) 

g + f 

-9 

11. 47 

-1. 500 

2.21 

+0 

6. 15 

-8. 5 

10. 83 

-1.398 

2. 12 

+0. 05 

6. 30 

-8 

10.20 

-1. 398 

2. 13 

+0. 10 

6.43 

-7. 5 

9. 56 

-1.204 

1.99 

+0. 15 

6. 55 

-7 

8. 93 

-1. 204 

2.02 

+0.20 

6. 65 

-6. 5 

8.29 

-1. 097 

1. 94 

+0.25 

6. 75 

-6 

7. 66 . 

-0. 954 

1.83 

+0. 30 

6. 85 

-5. 5 

7.02 

-0. 954 

1. 91 

+0. 45 

7. 02 

-5 

6.39 

-0. 900 

1.88 

+0.60 

7. 04 

-4. 5 

5. 75 

-0.750 

1. 80 

+0. 81 

6. 87 

-4 

5. 12 

-0.602 

1.72 

+1.10 

6. 31 

-3. 5 

4.48 

-0.602 

2. 00 

+ 1. 50 

5. 14 

-3 

3. 85 

-0. 500 

1. 97 

+1. 84 

4. 05 

-2. 5 

3. 23 

-0. 330 

1.91 



-2 

2.66 

-0 

1.81 
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